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Abstract. We prove the existence of a weak solution to a two-dimensional resonant 3x3 system of 
conservation laws with BV initial data. Due to possible resonance (coinciding eigenvalues), spatial BV 
estimates are in general not available. Instead, we use an entropy dissipation bound combined with 
the time translation invariance property of the system to prove existence based on a two-dimensional 
compensated compactness argument adapted from |37| . Existence is proved under the assumption that 
the flux functions in the two directions are linearly independent. 



1. Introduction 

This paper studies certain two-dimensional resonant 3x3 systems of conservation laws of the form 

k t = 0, l t = 0, 
ut + f(k,u) x + g(l,u) y = 0, 
which are augmented with L°° n BV initial data 

(1.2) k\t=o = k(x,y), l\ t= o = l(x,y), u| t=0 = u (x,y). 

The goal is to prove that there exists a weak solution to (|l.l|) - (jl.2|) . 
In recent years the one-dimensional version of the above system, 

kt = 0, 

(L3) ut + /(A,«)x = 0, 

has received a considerable amount of attention. This system may be viewed as an alternative way of 
writing a scalar conservation law with a discontinuous flux, namely 

(1.4) u t + f(k(x),u) x = 0. 

Equations like (|1.4p occur in a variety of applications, including flow in porous media, sedimentation 
processes, traffic flow, radar shape-from-shading problems, blood flow, and gas flow in a variable duct. 

If k(x) is a smooth function, Kruzkov's theory |22j tells us that there exists a unique entropy solution 
to the initial value problem for (|1.4[) , for general flux functions /. The scalar Kruzkov theory does not 
apply when k(x) is discontinuous. Instead it proves useful to rewrite (|1.4p as a 2 x 2 system of equations 
(11. 3p . which makes it possible to apply ideas from the theory of systems of conservation laws. 

As a starting point, it is necessary to introduce conditions on the flux f(k,u) that guarantee that 
solutions stay uniformly bounded. For example, one can require f{k,a) — f(k,b) — for all k, which in 
fact implies that the interval [a, b] C K becomes an invariant region. The system (|1 .4|) has two eigenvalues, 
namely Ai = and A2 = f u (k,u). Consequently, if f u (k,u) vanishes for some value of [k, u), then (|1.4[) 
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is nonstrictly hyperbolic and experiences so-called nonlinear resonant behavior, which implies that wave 
interactions are more complicated than in strictly hyperbolic systems. As a matter of fact, one cannot 
expect to bound the total variation of the conserved quantities directly, but only when measured under 
a certain singular mapping. A singular mapping that is relevant for (II. 3|) is 

*(fc,u) = J \f u (k,0\ d£. 

If {u p } p>0 is a sequence of "reasonable" approximate solutions of (jl.3|) . then one proves that the total 
variation of the transformed quantity z p :— ^(k,u p ) is bounded independently of p. Helly's theorem 
then gives convergence (along a subsequence) of z p as p J. 0. Since the continuous mapping u <— > ^(A^u) 
is one-to-one, u p also converges. 

A singular mapping was used first by Temple [40] to establish convergence of the Glimm scheme (and 
thereby the existence of a weak solution) for a 2 x 2 resonant system of conservation laws modeling 
the displacement of oil in a reservoir by water and polymer, which is now known to be equivalent to a 
conservation law with a discontinuous coefficient (see, e.g., [21] )■ Since then the singular mapping ap- 
proach has been used and adapted by great many authors to prove existence of weak solutions to resonant 
systems of conservation laws/scalar conservation laws with discontinuous flux functions, by establishing 
convergence of various approximations schemes (Glimm and Godunov schemes, front tracking, upwind 
and central type schemes, vanishing viscosity /smoothing method, . . .), see (the list is far from being 
complete) [BillGIlEJiaiaillliniin Similar ideas have been used 

also in the context of degenerate parabolic equations [16]. Regarding uniqueness and entropy conditions 
for scalar conservation laws with discontinuous coefficients, see jTTJ [TB] and the references therein. 

As an alternative to the singular mapping approach, the papers [TU [131 HE] has suggested to use the 
compensated compactness method and "scalar entropies" for the convergence analysis of approximate 
solutions. The results obtained with this approach are more general (and to some extent the proofs are 
easier) than those obtained with the singular mapping approach. 

All the papers up to now have addressed the one-dimensional case. The aim of the present paper is 
to take a first look at the multi-dimensional case, which is completely unexplored. More precisely, we 
will prove the existence of at least one weak solution to the initial value problem for the two-dimensional 
system 

Our existence proof is based on studying the " (e, 5) [ (0, 0) limit" of classical solutions u £ ' S of the 
uniformly parabolic equation 

uf + f(k s , u s >% + g(l 5 , u £ '% = e (uH + u^) , e > 0, 5 > 0, 

where k s , I s converge to k, I in L 1 1 oc (R 2 ), respectively, as 5 J. 0. 

Observe that we are essentially considering a scalar approximation scheme for see p| IT] ITC ] IT4" ] 

I15[ I18[ I41[ 142] for other scalar approximation schemes for one-dimensional discontinuous flux problems. 

Although spatial BV bounds are out of reach, we still have a time translation invariance property 
at our disposal, which, together with the assumption of BV initial data, implies that u t ' S is uniformly 
bounded in L 1 . Consider three functions F(k,u), G(l,u), H(k,l,u) defined by 

F*t = (fu) 2 , G u = (gu) 2 , H u = f u g u . 

We prove, at least under the assumption that e and S are of comparable size, that the two sequences 

F(k(x, y), + H(k(x, y),l{x, y),u £ > s ) y 

and 

H{k(x, y),l(x, y), u s >% + G(k(x, y), 
are compact in Wj~ c 1,2 (M 2 ), for each fixed t > 0. 
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The crux of the convergence analysis is then to prove that the above W,~ ' (R 2 ) compactness is 
sufficient to establish a "two-dimensional" compensated compactness argument in the spirit of the classical 
Tartar-Murat results for one-dimensional conservation laws [27j [28j [29j [38J [39] (see also [4]). Here we 
follow the recent two-dimensional compensated compactness framework developed in Tadmor et. al. |37j 
for nonlinear conservation laws. We extend their results to the case involving additional discontinuous 
"variable coefficients". Accordingly, we make the nonlinearity assumption that for each fixed k,l the 
functions u t— > f u (k,u) and u i— » g u {l,u) are almost everywhere linearly independent (see (|2.4p in the 
next section for a precise statement). Our main existence result is based on an application of the two- 
dimensional compensated compactness lemma with "variable coefficients" — lemma [3~2l stated in Section 
Obelow. Granted the nonlinearity assumption, it then yields that (a subsequence of) u e ' S (•, •, t) converges 
in L 1 1 oc (R 2 ) to a bounded function u(- } -,t), for a.e. t > 0. Since u £,s is uniformly L 1 Lipschitz continuous 
in time we obtain, in Section 2] below, our main Theorem 12. f 1 stating that u £ ' 5 — > u in L 1 1 oc (R 2 x R + ) 
and that the limit function u is a weak solution of Ql.l|) - (|1.2|) . 

Although we have chosen to analyze the vanishing viscosity/smoothing method, the techniques used 
here for that purpose can also be applied to various numerical schemes, including appropriate two- 
dimensional versions of the scalar finite difference schemes studied in ITU HH1 EH H21 . 



(2.2) 



(2.3) 



2. Assumptions and statement of main results 

We start by listing the assumptions on the initial conditions uq and the fluxes k,l, /, g that are needed 
for the existence result. 

Regarding the initial function we assume 

(2.1) u a G L°°(l 2 )nW(l 2 ), a<u <b for a.e. in R 2 . 

For the discontinuous coefficients k, I : R 2 — > R we assume 

jk,le L°°(R 2 )n W(R 2 ), 
\a<k,l</3 a.e. in R 2 . 

For the flux functions /, g : [a, 0\ x [a, b] -^Iwe assume 

{u i— > f(k,u),u i— * g(l,u) G C 2 [a,b] for all k,l G [a,/?]; 
k i— ► f(k,u),l i— > g(l,u) G C^a,/?] for all u G [a, b]. 

Moreover, we make the nonlinearity assumption which excludes the possibility of £if(k,u) + £,2g(l,u) 
being an affine function (in u) on any nontrivial interval for all k, I G [a, /?], 

V |£| = 1 and k, I G [a, 0\ : £if(k, •) + ^2g(l, •) M affine function on any nontrivial interval. 

In its slightly stronger version, this assumption requires that f u (k, •) and <?«(£, •) are a.e. linearly inde- 
pendent so that the symbol s(£, k, I, u) :— £i/«(fc, u) + £,2gu(l, u) satisfies 

V |^"| = 1 and k,l€[a,f3]: meas{w | s(£, k, I, u) — 0} = 0. 

This is a straightforward generalization of the notion of nonlinearity found in |25j , in their study of kinetic 
formulations for nonlinear conservation laws. 

Finally, we need to know that our approximate solutions stay uniformly bounded. For example, this 
is ensured by the assumption 

(2.5) f(k, a), f(k, b), g(l, a), g(l, b) = for all k, I G [a, 0], 

which implies that the interval [a, b] becomes an invariant region. Of course, one can relax assumption 
(|2.5[) . A sufficient condition for the invariance of the interval [a, b] is that the divergence of the vector 
field (x,y) i—* (f(k(x,y),u),g(l(x,y),u)) is nonnegative when u — b and nonpositive when u — a. Let 
us emphasize that an assumption like (|2.5p is essential to our analysis; Without it solutions can possess 
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concentration effects, which is a weil-known feature of, for example, linear transport equations with 
discontinuous coefficients. 

We are now ready to state our main result. 

Theorem 2.1. Suppose (|2.1[) . (|2.2p . (12. 3|) . (12. 4p . and (|2.5p /ioZ<£ Then, there exists a weak solution of 
the initial value problem (fTTj) - (TP)) . u e L°°(R 2 x R + ) n Lip(R + ; L 1 (R 2 )), satisfying 



(u4> t + f(k(x, y), u)4> x + g(l(x, y), u)4> y ^j dx dt 
+ / u {x)(j){x,0)dx = 0, € 2?(R 2 x [0, oo)). 



77ie weafc solution, u, can be constructed as a strong L 1 1 oc (R 2 x R + ) -limit of classical solutions u E,s of 
uniformly parabolic problems, 

(2.6) uf + f(k s ,u s > s ) x + g(l 5 , u E ' 5 ) y = eAu"' 5 , 

with the smoothly mollified coefficients, k s :— ujs * k and I s :— tog * I (outlined in section \2.1\ below). 

The proof of this theorem is given in the following two sections. Remark that the BV assumption on 
the coefficients k, I made in (|2.2|) is used twice in this paper. First, it is used to prove Lipschitz regularity 



in time, in lemma l4~2l below: then, we use it to prove Wi oc ' (R 2 ) compactness of the entropy production 



loc 

for each fixed t > in lemma 14.41 below. 

We close this section with the following summary. 

Remark 2.2. Stated differently, Theorem 12.11 shows that there exists a weak solution to the following 
two-dimensional scalar conservation law with discontinuous coefficients k, I £ L°°(R 2 ) n BV(M. 2 ): 

u t + f{k(x, y), u) x + g(l(x, y), u) y = 0, 

«|t=o = u a e L°°(R 2 ) n BV(R 2 ). 

We note in passing that the solution operator in this case of discontinuous "variable coefficients" is not 
translation invariant in space and hence the ^-contraction property of (|2.6p does not imply spatial BV 
compactness. 

Moreover, if we let uq(-) >— > u(t, •) denote the mapping of (|l.ip - (|1.2p . so that u(t,-) is a (vanishing 
viscosity) weak solution constructed in theorem [23] then by adapting standard arguments, we can prove 
that the mapping is compact with respect to the L\ norm. 

3. A COMPENSATED COMPACTNESS LEMMA 

In this section we prove a " two-dimensional" compensated compactness lemma. We refer [8j [26l [27l 
|2"51 [2^1 13"51 [55] for background information on the compensated compactness theory. We start by recalling 
the celebrated div-curl lemma. 



Lemma 3.1 (div-curl). Let fl C R 2 be an open domain and let p > denote a parameter taking its 

}p>0> 



values in a sequence which tends to zero. Suppose DP -± D, EP -± E in (L 2 (n)) and {dbr DP} ^ 
{curli? p } p>0 belong to a compact subset ofW^ c ' (Q). Then, after extracting a subsequence if necessary, 
we have D p ■ E p -^~D - ~E in D'(fi) as p \ 0. 

The compensated compactness lemma below is tailored for two-dimensional equations, whose spatial 
part involve discontinuous coefficients: 

f(k(x, y),v(x, y)) x + g(l(x, y), v(x, y)) y . 
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If g(l,u) — g(u) and f(k,u) — f(u), then the lemma below coincides with the two-dimensional result 
of [23 Theorem 3.1]. If we set g = then the result coincides with Tartar's compensated compactness 
lemma for the one-dimensional scalar conservation with genuinely nonlinear flux /. 

Lemma 3.2 (Compensated compactness). Let fl c K 2 be an open domain. Let k,l,f,g be functions 
satisfying (|2.2| . (|2.3| . (|2.4| . and (|2.5p . Suppose {v p (x, y)} p>0 * s a sequence of measurable functions on 
Q that satisfies the following two conditions: 

(1) There exist two finite constants a < b independent of p such that 

a < v p (x, y) < b for a.e. (x, y) G fi. 

(2) Let the functions F, G, H be defined by 

F u (k,u) = (f u (k,u)) 2 , G u (l,u) = (g u (l,u)) 2 , H u (k,l,u) = f u (k,u)g u (l,u). 
We assume that the two sequences 

{F (k(x, y),v") x + H (k(x, y),l(x, y), v")\ , 
(3.1) J ">° 

[H (k(x, y),l(x, y),v") x + G (l(x, y), v p ) y \ p>Q 

belong to a compact subset of ' (1 x 
Then, there exists a subsequence of {v p (x, y)} p>0 that converges a.e. to a function v £ L^^M 2 ), and 
a < v(x, y) < b for a.e. (x, y) € CI. 

Proof. To simplify the notation let 

FP :=F(k(x,y),vP), G p := G (l(x,y),v p ) , H p := H (k(x, y),l{x, y),v p ) , 
and denote their L°°(R x R + ) weak-* limits by F, G, H, respectively. Introduce the vector fields 

D p =(f p 1 H p ], E p =(-G p ,H p ^ 



and denote their respective L°°(IR x R + ) weak-* limits by D, E. 

Thanks to ()3.1|) . we can apply the div-curl lemma to the sequences {D p } p>0 , {E p } p>0 to produce 



D-E = D-E a.e. in ft, 

that is, 



H 2 -FG=(H) -FG, 

which implies 



(3.2) (H-H)-(F-F)(G-G)=0. 

Fix c = c(x, y) E L°°(£l). Following [37], we now consider the function / : [a, b] —> M. defined by 

I(v) = (h (k(x, y),l(x, y),v) - H (k(x, y),l(x, y),c)^ 

- (f (k(x, y),v)-F (k(x, y), c) ) • (g (l(x, y),v)-G (l(x, y), c) ) . 

Note that 

I(v p )= {[H p -H] + [H-H(k(x,y)J(x,y) 7 c)]y 

-([F p -F] + [F-F(k(x,y),c)]) ■ ([ G p - G] + [ G - G (l(x, y), c) ] ) . 
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Using this and (|3.2p , we compute 
(3.3) = (H - H (k(x,y)J(x,y),c)y - (T - F (k(x,y),c)y(G ~ G(l(x,y) lC ) 

By the Cauchy-Schwarz inequality we have for any u 6 [a, b] 

(h (k{x, y),l(x, y),u) — H (k(x, y),l(x, y),c)j 

( 3 - 4 ) = fu(k(x, y), 09u(l(x, y), d{\ 

< (F (k(x, y),u)-F (k(x, y), c) ) (g (l(x, y),u) - G (l(x, y), c)) , 

and hence /(•) < with 1(c) = 0. Thanks to (12. 4p . the Cauchy-Schwarz inequality in (|3.4I) is in fact a 
strict inequality. This shows that the function I(v) has a strict global maximum at v = c. 
Since u <—> F (k(x, y), u) is strictly increasing, we can choose c as 

c(x,y) := F^ 1 (k(x,y), F(x,y)) , 

so that (|3.3[) becomes 

H^) = (H-H(k(x,y),l(x,y),c)y. 

Since /(•) < 0, we conclude that H = H (k(x,y),l(x,y),c), and thus I(v) = 0. In fact, we have 
I(v p ) -> a.e. in Cl. 

Using the fact that c is a strict maximizer of I(v), we have 

I( v ) — ~C a , whenever \v — c\ > a, 

for some constant C a > that depends on a. Consequently, 

meas {\v p — u\ > a} < — - 11 I(v p (x,y)) dx dy — ► as p J, 0. 

Got J Jfln\vP-c\>a 

Since a > was arbitrary, v p — > c in measure, which in turn implies that a subsequence of {v p } s>0 
converges to c a.e. in Cl. □ 



We remark that the idea of using the Cauchy-Schwarz inequality along the lines of (13. 4p for proving 
strong compactness can be traced back to [HH 03] 

4. Proof of Theorem 12. II 

Let k s ,1 s , Ug be smooth functions converging strongly to k,l,uo respectively. More precisely, let us € 
Cq°(R) be a nonnegative function satisfying 



uj(x) = for |x| > 1, / / cj(x) dx = 1. 

J Jm. 2 

For 5 > 0, let uu$(x) — iu; (|) and introduce the mollified functions 

k s =ujs*k, 1 S — lus*1. 
We approximate the initial data uq by cut-off and mollification as follows: 

where x$( x ) = 1 f° r \ x \ ^ an< ^ otherwise. In particular, we have the estimate 

+ K) J L1(R2) < J// r2 (N)J + d ^ ^ ] i w °iw (K2) 
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Observe that for h s = k s , I s , u s and h = k,l, u , we have h s € C°°(R 2 ) and 

h s -> h a.e. in M 2 and in L P (R 2 ) for any p € [1, oo) as S I 0. 

Additionally, Uq is compactly supported. 

Now let u e ' 5 be the unique classical solution of the uniformly parabolic equation 

(4-2) uf + f(k 5 ,u £ > 5 ) x + g(l s , u £ ' S ) y = e + vtf) , 

with initial data u £ ' S \t=o = Uq. The proof proceeds through a series of lemmas, which in the end show 
that for each t € [0, T] a subsequence of u £ ' S (-, •, t) converges a.e. as e, <5 | 0. 
Our first lemma confirms the uniform bound. 



Lemma 4.1 (L°° bound). There exists a constant C > 0, independent of e,5, such that 

\l°°(r 2 ) 



\u^(;;t)\\ Loo(K2 .<C, for all t€(0,T). 



Proof. The proof is standard and exploits assumption (|2.5[) to conclude that a < u £ -' 5 (x,y,t) < b for 
a.e. (x,y,t) eR 2 x R + . □ 

Using that (|4.2[) is translation invariant in time, we can prove that u e t ,S is uniformly bounded in 
L°°(M.+ ;L 1 {R 2 )). 

Lemma 4.2 (Lipschitz regularity in time). Suppose the two smoothing parameters e and 5 are kept in 
balance in the sense that 

(4.3) 5 = Ce, for some constant C > 0. 

There is a constant Co (which is possibly dependent on uq but otherwise is independent ofe,5), such that 
for any t > 

ff \d t U £ > S (;;t)\dxdy<C . 

Proof. To prove this, set w £ ' S — u £ t ' S . Then w satisfies 

wf + (f u (k s ,u £ ' s )w £ >% + (g u (l s lU £ - s )w £ - s ) x = e {w £ J + w £ yy s ) . 
Multiplying by sign(w e '' 5 ) gives, in the sense of distributions, 

\w £ ' s l + (f u (k s ,u £ - 5 ) \w £ ' s \) x + (g u (l 5 ,u £ ' s ) \w £ < 5 \) y 



e (\w £ ' s \ + \ w £ - s \ I - esigr 



5nV tf )(K' a ) 2 +«) 2 ), 
since f u (k s ,u £ - 5 )w £ - s sign' {w £ - s )w £ x 5 = g u (k s ,u £ ' S )w £ ' S sign' (w £ ' S )w £ y ' S = 0. Hence 

j t JJ^\w £ - s \(x,y r )dxdy <0, 
which, due to ()4.1j) and (|4.3j) . concludes the proof. □ 
Thanks to the previous lemma, we also have uniform L°°(R + ; L 2 (M. 2 )) control over ^feu 6 :" and y/euy . 
Lemma 4.3 (Entropy dissipation bound). There is a constant C, independent of e,5, such that 

((u £ x S (;;t)) 2 + (u £ y s (-,;t)) 2 ) dx dy < C, for any t > 0. 
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Proof. Multiplying (|4.2[) by u £ ' 5 and then integrating yield 
r(«' 5 (-,-,t)) 2 + (^(.,.,i)) 2 )dxdy 



uu e t ' S dx dy 



h(k s ,Odnk s x (x,y) 



dx dy 



dx dy. 



In view of Lemmas 14.11 and 14.21 and the BV regularity of the coefficients, we derive easily the uniform 
bound 



< C 



({uf(;;t)f + (uf (;;t)f)dxdy 

+ Ww(R 2 ) + \1\bv(r 2 )^J ' * > 0, 





e,<5 







L°°(R + ;L 1 (M 2 )) 
e,6\ 



for some constant C that is dependent on ||w e ' ||_l~( K 2) but otherwise is independent of e,S. 



□ 



Lemma 4.4 (Pre-compactness at each time instant). Suppose the two parameters e and S are kept in 
balance in the sense that (|4.3p holds. With F , G, and H defined in Lemma \3.2l the two sequences 



{F (k(x, y),u^ s ) x + H (k(x, y),l(x, y), u e > s ) y ) ^ ^ 

[H (k(x, y), l(x, y), u s < 5 ) x + G (l(x, y), u e ' a )„} e g>Q 

then belong to a compact subset of WJ7 ' (K 2 ), for each fixed t > 0. 

Proof. Let <\> = (j>(%, y) € 2?(R 2 ), and, for each fixed t > 0, introduce the distribution 

4>) = [[ (F (k(x, y),u s ' s ) cf> x + H (k(x, y), l(x, y)u^ s ) y ) dx dy. 



Let us first write C e ' S — C e { S + C^f , where 



e,S 



e,S 



(F (k(x, y),u^ s ) - F (k 5 (x, y), z/' 5 )) <j> x dx dy 
+ II (H (k(x, y), l(x, y), u^ s ) - H (k\x, y), I s (x, y),u^))c^ y dx dy, 
(F (k 5 {x, y), u'- s ) <p x + H (k 5 {x 7 y),l 5 {x 7 y), u s < 5 ) <j> y } dx dy. 



C 



In what follows, we let fl denote an arbitrary but fixed bounded open subset of R 2 . With <f> S W ' (Q), 
we have by Holder's inequality 



L 'l ) 



<C{\\k-k 6 



\k~k 6 



\L 2 {n) 1 II llL 2 (n) 
as S i 0. Thus, •! C e { 5 \ is compact in W~ 1,2 {VL), for each fixed t. 

L J £.<5>0 



lw l,a (n) — * °> 
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Multiplying 



by f u (k s (x,y),u e > 5 ) yields 



F (k 5 , u^ 5 ) x + H {k 5 , I s , u^ 5 ) v = if + if + If + If + I, 



re, 5 
l 5 > 



where 



f = (suff u (k s ,u s > s )) x + (euff u (k s ,u s ' s )) y , 

f = -e (uf) 2 f uu (k 5 ,u^ s ) - e (uf) 2 f uu (k 5 , u^ s ), 
f = -euf f uk (k s ,u s ' s )k s x - euff uk (k 5 ,u^)k s y , 
f = F k (k s , u e < 5 )k 5 x + H k (k s ,l 6 , u"' s )k s y + Hi(k 5 , I s , u 6 ' 5 )l s y 
" -- -f u (k^u e ' 5 )f k (k^u s < 5 )ki-f u (k^u s ' s ) gi (l\u s ' s )l" 



Hence, there is a natural decomposition of Cf into six parts. We name the six parts Cf , ^221 

•'-'2,3' ^2,4' ancL J "2,5- 



Regarding Cf , 



C 



< C 



e,S 



H0ll£-(n). ^eCo(ft), 



which yields 



C 



< C, where M(Q) = ((7 C (0)) denotes the space of bounded measures on Q 



M(Q) 

In view of Lemma 14.31 and the uniform boundedness of the solutions, it follows that 



e,5>0 



compact (and in fact converges to zero) in W 1,2 (fi) and 
Next, for any (j> £ C C (U), observe that 



C 



< C, for each fixed t > 0. 



lf(j) dx dy 



< C 



£ [VZ 



5 ) dx dy 



C 



\ek y \ \k y \ 



lf<t> dx dy 



e (ufY dxdy j 

The point here is to have e and S in balance, so that we can ensure \ek 5 x | , \sk y \ < C. More precisely, we 
have \ek s x | , \sk y \ < C|, and by choosing e,S according to (|4.3p we achieve this balance. Consequently, 

< C, for each fixed t > 0. 



and thus 



-2,3 



Finally, using the BV regularity of the coefficients and the boundedness of the solutions, 

<C\ 



lf(j) dx dy 



lL°°(fi) 



and thus 



C 



Similarly, 



< C, for each fixed t > 0. 

M(n) 

Il^lf^dxdy < C\\(f)\\ Lao , Q y and thus 



C 



M(n) 



< C, for each fixed t > 0. 



Summarizing, we have shown that the sequence of distributions {£ £ ' S } e S>Q satisfies the following two 
properties: {i} each distribution is the sum of two terms — one is compact in W _1,2 (0) and the other 
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one is bounded in M.(fl). In addition, Lemma HTTl implies that {ii} {£ e,<5 } e g belongs to a bounded 

subset of VK~ 1 '°°(ri). We now appeal to Murat lemma [2§]> which guarantees that {C £ ' S } £ S>Q belongs to 

a compact subset of W _1 ' 2 (f2). This concludes the proof of the first part of the lemma, since 51 was an 
arbitrary bounded open subset of R 2 . The second part of the lemma can be proved in a similar way. □ 

Concluding the proof of Theorem \2.1\ By combining Lemmas 13.21 and 14.41 we conclude that u 6 ' S (-, •, t) is 
pre-compact a.e. for each t £ [0, T]. Together with a diagonal argument, we can prove that u £ > s (-, -, t) 
converges along a subsequence a.e. in R 2 and in L 1 1 oc (R 2 ), for each fixed t > 0. Lemma l4~2l implies that 

||u £ ' 5 (-,-,< + r)-^(.,.,i)|| il(R2) <CV, Vre (0,T-t), 

and using this L 1 time continuity estimate it takes a standard density argument to show that there exists 
a subsequence of {w e,<5 } g S>Q that converges to a limit function u a.e. in R 2 x R + and in L 1 1 oc (R 2 x R+). 

Moreover, the limit u belongs to ^(R 2 x R + ) n Lip(E + ;£ 1 (R)). 

Equipped with the strong convergence it is easy to prove that the limit u is a weak solution. This 
concludes the proof of Theorem 12.11 □ 

We close the paper by a couple of remarks. First, we note that the apriori bounds in Lemmas 14.1) 14.21 
and !4.31 being uniform in time, enabled us to deduce pre-compactness at each fix t > 0, thus circumventing 
the temporal argument required in [371 Appendix A]. Second, we have herein exclusively dealt with 
problems that are spatially two-dimensional. A possible strategy for going beyond two dimensions is to 
apply the compactness framework of Panov pH] [S3 H2] • 
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